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1 Introduction 



Heavy-light bound states in QCD can be described efficiently by an expansion in the 

inverse heavy quark mass. Already in the early days of the associated effective field 
theory, HQET [1-3], the mass splitting between vector and pseudo-scalar heavy- light 
mesons served as a phenomenological argument for the absence of large higher order 
corrections in the expansion. 

Consider QCD with A'f light quark flavors and a heavy flavor, the b-quark, iVf = 4 
being the case realized in Nature. The splitting 

Am^ = m|* — (1-1) 

then has an asymptotic behavior for large quark mass mb which is characterized by one 
renormalization group invariant (RGI) observable, 

4Xf^' = hm |[26o5'(mb)]-"°/''°Amn, (1.2) 
(7o = 3/(87r2), 6o = (11 _ 2^,)/(i67r2)) 

of dimension [mass]^. Since the limit exists, this quantity is uniquely defined in QCD. 
In the above the definition of is irrelevant as long as it is renormalized at a scale of 
order mb- 

As a rather non-trivial statement, the effective field theory predicts that A2^*^^ can 
unambiguously be computed in that theory, where the b-quark is treated as static. 
There it is expressed as an expectation value 

A^i = l{B\Of^^\B) / {B\B) , (1.3) 
O^^ = limj2bo-9\fi)r''^''''' O^^M , (1.4) 

of in the zero-momentum B-meson state \B). The operator 0^;^ is related to the 
bare local operator 

Cspin(a;) = iph{x)^.Fki{x)akiiph{x) = iph{x)o- ■B{x)ip^{x) (1.5) 

by a multiplicative renormalization depending on the adopted scheme S and a renor- 
malization scale fi - but neither depends on a scheme nor on a scale. 

The matrix element of the bare operator can be computed non-perturbatively by 
lattice simulations of HQET [4-6]. As stated in these references, a significant source of 
uncertainty remained in the connection between the bare operator and the RGI one (or 
one renormalized in the MS scheme), which has only been established perturbatively [7, 
8]. This uncertainty made it impossible to decide whether the splitting is significantly 
underestimated in the quenched approximation (iVf = 0) or not. 

In this paper we develop the non-perturbative renormalization of Ogpin- We follow 
the general strategy of the ALPHA-coUaboration [9-11], specialize it to the operator in 
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question and perform an explicit computation in the quenched approximation. In partic- 
ular we define a suitable scheme using Schrodinger functional (SF) boundary conditions, 
and compute the 2-loop anomalous dimension in this scheme. We then evaluate the scale 
dependence of the operator non-perturbatively for iVf = 0, between fi ~ 0.3 GeV and 
fi 100 GeV. Using the high energy end of the results and eq. (|1.4p supplemented 
with the 2-loop anomalous dimension, the connection to the RGI operator is realized. 
Finally the total Z-factor between bare and RGI operator is obtained for the Wilson 
gauge action and several HQET discretizations. Readers solely interested in the final 
result for the Z-factor may find it in Sect. 15. li 

For a comparison to the experimental mass splitting at finite mass, it is important 
to include radiative corrections beyond the 1-loop ones incorporated in eq. (II. 2|) . We 
do this in the form 

Am^ = 2^^^5_±!!^C'spin(Mb/AMs)A^^' + 0(l/mb) , (1.6) 

written in terms of RGI's with a function Cspin known up to corrections 0(a^(mb)) ~ 
4% and discussed in some detail in Sect. [21 Obviously the perturbative uncertainty 
can be estimated more reliably and reduced by a higher order continuum perturbative 
computation in QCD. We note that the final renormalization factor also applies 

to spin-dependent potentials computed in lattice gauge theory (see Sect. [3]). 

The reader is not to confuse the present approach with the one of [12,13], where, 
through a non-perturbative matching between QCD and HQET, also functions such 
as Cspin are determined non-perturbatively. While in general the strategy of [12, 13] is 
essential, the more traditional path is viable here because Ogpin does not mix with lower 
dimensional operators. 

Before entering the discussion of the renormalization of Ogpinj we briefiy address the 
question of the precision that can be expected from eq. (jl.6p . For this purpose we boldly 
also treat Ami = "^d* " "^D HQET. So A'f = 3 in all places. With a iVf = 3 QCD- 
parameter of A^^ = 300(100) MeV and with = 1.55 GeV, Mb = 6.69 GeV [14] we 
find Mc/Mb ^ 0.23 and Cspin(Mc/A^jg)/Cspin(Mb/Aj^) = 0.94 . HQET then relates 
the splittings as Am'^/Am^ = 1.41(2), where the uncertainty is due to the generous 
error in A^^. With the quenched input values Aj^ = 238 MeV, = 1.65 GeV, Mb = 
6.76 GeV [11,13,15,16] this ratio changes only slightly, namely to Am^/Am^ ~ 1.44. 
This is to be compared to Am^/ Am? = 1.14 from experiment. 

Since the charm mass is only moderately large, such a 25% deviation is not unex- 
pected. Scaling this correction to the B-system, we expect an accuracy of order 5-10% 
for the HQET prediction of Am^. Earlier quenched approximation estimates with per- 
turbative renormalization found values for Am^ which were lower than the experimental 
number by between 50% [4,5] and 20% [6]. Renormalizing the same matrix elements 
non-perturbatively we will find the difference to experiment significantly reduced in 
Sect. El 
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2 HQET and A^^^ 



2.1 Lattice action 

We briefly define the effective theory in a lattice regularization, using the notation 
of [17, 18]. The heavy quark fields are taken to have 4 components with the constraint 

P+V'h = V'h, V^h^+=V^h ^'+ = ^(1 + 70). (2.1) 
With the lattice backward derivative 

D^Mx) = I [M^) - W\x - ad, 0)Vh(x - aO)] , (2.2) 

and the mass counterterm dm-^j^, the static action (i.e. lowest order HQET) is written 
as 

= E V^h(x)(^^ + 5mw)V'h(x) . (2.3) 

1 + a omw ^ 

Different gauge connections W have been found to be very useful to improve the statisti- 
cal precision in numerical simulations [18,19]. They play a role only when we discuss the 
non-perturbative results. Until then the reader may think of W{x,Q) as the standard 
timelike link. In fact that choice defines the original Eichten-Hill action [1]. 



2.2 Conversion functions 

The operator C'spin(x) in eq. (jl.Sp is given in terms of the flelds entering eq. (j2.3p with 
the normalization specified there. The lattice version F^y of the gauge field tensor is 
defined by the clover leaf term, see e.g. [20]. Ogpin appears as a first order correction 
in 1/mb in HQET and induces the spin splitting. Usually, the splitting is written in a 
form different from eq. (jl.6p . We want to briefly explain why we choose the latter. 
The more common form is 

ml. -ml = 4C;jJ,^*^^mb)A2(mb) + 0(l/mb) (2.4) 
A2(mb) = \{B\dS,{^i = m^)\B) (2.5) 

where and mb are renormalized in the MS scheme. This is arrived at by starting 

from the formal expression 

mB* -me ~ |^(S|Ospm|5) / • (2.6) 

One renormalizes Ospin in the MS scheme, identifies the mass mb with the (perturbative) 
pole mass mQ,b and defines the remaining factor as a matching coefficient C'™|'g^^(mb). 
Finally one uses 2mQ b = m-Q* + m-Q + 0(A), dropping the 0(A) correction. 

The matching coefficient C™^g"^'^(mb) = 1 + Ci(f'{m\^) + ... is independent of the 
particular matrix element. S 4GeV is reasonably large and there is no mixing 
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Figure 1: Conversion functions for A'^f = 0. Dashed lines for 1-loop anomalous dimension and contin- 
uous line for 2-loop 7. In the latter case the parametric uncertainty is 0(a^(mb)). The abscissas of the 
c quark and the b quark are marked by vertical dots. 



with lower dimensional operators, C^^^^ can be approximated by perturbation theory. 
It is known including the Ci^^(mb) term [7]. 

In the above form, the matrix element and the HQET parameter A2 depend on the 
arbitrary renormalization scheme (MS). Such a spurious dependence is easily removed 
by introducing RGIs (see e.g. [21] and Sect.III.3.1.of [22]), in particular O^^, eq. (fLll) . 
It is related to the bare operator Ogpin in a particular lattice regularization via 

C^?n = ^sph?(9o)Ospin • (2.7) 



We now have 



ml. -ml = 4C„,ag(Mb/Aj^)Af^i + 0{l/m^) 



with a function Cniag(-^b/A^) written in terms of the RGI mass of the b-quark and 
the QCD A-parameter. Using existing perturbative computations [7, 23-25] it is easily 
evaluated by integration of the RG equations (see e.g. [21]). 

The comparison of the successive perturbative approximations for Cmag in Fig- H 
indicates a relatively large perturbative error even at the 2-loop level (for the anoma- 
lous dimension 7). At a scale of 4GeV this is somewhat unusual. However, one can 
understand this behavior by noting that the definition of the matching factor C^^^^^ 
involves the pole quark mass, which is unphysical. That mass has a perturbative rela- 
tion to short distance masses such as rnjj^ which is badly behaved (large perturbative 
coefficients). Also the relation pole mass to RGI mass has this property. 

Our eq. (jl.6p follows from choosing directly the RGI mass instead of the pole mass 
in eq. (|2.6p . In perturbation theory one then has the relation 

Cspin(Mb/A5^) = i^C^,g(Mb/Aj^). (2.9) 
mQ,h 
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The function Cgpin is reproduced from [21] in Fig. [TJ One notices that the two available 
perturbative approximations are much closer than they are for Cmag- We expect this to 
hold also at higher orders of perturbation theory. Still, since the very close agreement of 
the two approximations for Cgpin may be accidental, we will use an error of a'^{m]:,) ss4% 
at the mass of the b-quark. This perturbative error is currently being reduced by an 
explicit 3-loop computation [26]. 

In summary, the form eq. (jl.6p is written in terms of properly defined RGI's and 
Cspin appears to have a well behaved perturbative expansion. This is a good basis for a 
computation of the mass splitting. 

3 Schrodinger functional renormalization scheme 
3.1 Definition 

We want to formulate a renormalization condition for Ospin in a finite volume, which 
allows us to carry out a non-perturbative computation of the associated renormalization 
factor Z^^, following the general strategy of [11]. We choose Schrodinger functional 
boundary conditions since this allows us to perform accurate simulations and also per- 
turbative computations; see [22] for a review. There is an additional reason for this 
choice. With any kind of periodic boundary conditions, any correlation function with 
Cspin vanishes at tree-level. In order to avoid this, we choose boundary conditions which 
induce a non-trivial background field F^i, / at tree level. This ensures a good sig- 
nal in the MC simulations at weak coupling and means that a 1-loop computation is 
sufficient to compute the Z-factor up to and including 0{g'^). Since the operator does 
not contain any light fermion fields, we avoid these altogether in the definition of the 
correlation functions. For A^f = we then end up with a pure gauge theory definition 
(without valence quarks). Furthermore one easily sees that the possible dimension six 
operators which are necessary for the 0(a)-improvement of Ogpin do not contribute here. 
As is well-known, there is also no mixing with Okin = V'h^^V'h) the other dimension-5 
operator of hqetQ 

These considerations motivate the following choice. We take an Lq x Li x L2 x L3 
geometry and adopt Dirichlet boundary conditions inducing a background field as in [27]. 
But we choose Dirichlet conditions in the 3-direction, 

C/(x, /x) 1x3=0 = exp(aC) , C/(x, /x)|x3=L3 = exp(aC') , /i = 0,1,2, (3.1) 

keeping periodic boundary conditions with respect to xo,xi,X2 (remember that time is 

^The static action is invariant under the space-dependent transformation 5ipii{x) — 
a;(x) ^h{x) , 5tf}^^{x) = — cij(x) tpi^{x) , which corresponds to the conservation of the local b-quark num- 
ber [17]. While Ospin is invariant under this symmetry transformation, Okin is not. 
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already distinguished in the static quark action). The Abehan fields 



1 1 
C = -diag((/)i,(?;)2,(A3) = ydiag(-7r/3,0,7r/3) , 

^ ^ (3.2) 

C = -^diag((/.;,(/>'2,(/.3) = -^diag(-^,7r/3,27r/3), 

of "point A" [27] are chosen and we set L = Lq = Li = L2. The classical solution then 
has non- vanishing F^s, independent of the space-time position. Note that the strength 
of the fields scales with L. 

A natural choice of a renormalization condition is then 



.SF ,,,^'(5l(x + |6)Ospin(x)) _ L^Si{x + ilQ)0,^Ux)\ 



Z^^ (L) — 
'P'"^ {Si{x+^0)Si{x)) {Si{x + ^6)Si{x 

with X3 = L3/2 = L/2. The spin operator 



(3.3) 



90=0 



^kix) = — i;^{x)akWHx - ad, 0)Vh(x - aO) (3.4) 

is the simplest choice to obtain a non-vanishing trace in spin space. It is the (local) 
Nother charge of the invariance of the action under the transformation 

Siph{x) = Lj(x)akTph{x) , Si/j^ix) = -u;(x)tp^{x)ak , cjfc = -^eijkcrij , (3.5) 

with infinitesimal space-dependent parameter u;(x) and [(7fc,c7;] = iekimf^m- The Ward 
identities derived from this invariance imply that Sk{x) is not renormalized. Thus no 
additional factors are necessary in eq. (j3.3p . 

We have formulated the renormalization condition in terms of correlation functions 
of local operators to make it obvious that the standard theory of renormalization in- 
cluding 0(a)-improvement applies. However, integrating out the static quark fields we 
arrive at a form which is more natural for explicit computations, perturbative and non- 
perturbative. This step also shows immediately the connection to the spin-dependent 
potentials. With the explicit form of the static propagator [18], one obtains 

(5i(x + fO)5i(x)) (Tr(Po(a;))) 
where the Polyakov loop operator 

Vf^ix) = W{x, ^i)W{x + afi,ij)... W{x + (L - a)[i, n) (3.7) 

enters. Here F23ix) stands for the clover leaf discretization of the field strength ten- 
sor [20] (an alternative discretization will also be considered in the non-perturbative 
computations). Now the renormalization condition eq. (j3.3p is given in terms of the 
expectation values of a (traced) Polyakov loop and of a (traced) Polyakov loop with the 
insertion of a S field. 
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At this point it is useful to digress for one paragraph in order to exhibit the relation 
to spin-dependent potentials [28,29]. In standard notatiorH they can be defined as 
(periodic boundary conditions in all directions, = + + 

'-ifv^ir.,) = lim „y^(Tr(n(0)B.(0))TV(n(x)tB,(.))> 

r2 3 ^ '^^ ^ spinwpyj ^^^^^ ^ (Tr-Po(0)Tr-Po(a;)t) 

xq 

Vr{r,^^)-Vi^{r,f^) _ .^SF ^ /,Al2 „ ( Tr(Po(0)i?i(0))TY(Po(x)t j?i(x))) 

Translating from the potential renormalized in the SF scheme at renormalization scale fi 
to an RGI potential works just as explained in Sect.[2l In contrast to the static potential 
where a difficult to determine additive renormalization results from draw-, there is no 
additive renormalization in the above equations; 6m\\r drops out. Note that we have 
nothing to add to the phenomenological relevance of such potentials [30]. Rather we 
remark that these objects, computed recently in [31], have to be renormalized with 
Zgpin, which arises in HQET. 

Returning to our renormalization condition, we note that it is natural to use the 
equivalence of all coordinates in Euclidean space to switch to the usual SF boundary 
conditions, where 



L^(Tr(P3(x)j^i(x))) ^ L2(TV(P3(^)i^i(x))) 



7SF 



' at xo = ^ (3.8) 

90=0 ^ 



El (x) = iFqi (x) , Dirichlet boundary conditions in time. 



From now on we retain these boundary conditions and this coordinate system. The 
tree-level value is 

L2(Tr(p3(x)Ei(x))) ^1 + ^3 



{TtVsix) 



+ 0{{a/Lf). (3.9) 
go=o o 2 - V3 



Corresponding formulae for finite a/L, which are used in the non-perturbative definition 
in order to assure Zspin = 1 at tree-level, are given in eqs. ()A.34l IA.35[) . 

3.2 The 2-Ioop anomalous dimension 

Our strategy for computing the RGI renormalization is to non-perturbatively evaluate 
■^spin(l//") up to /i = O(lOOGeV) and then evaluate O^^VC'spm(/u) in perturbation 
theory. In the latter step the anomalous dimension 7(g), defined by the renormalization 
group equation 



d_ 



/x^OSpF„ = 7(5SF)Ofp'^n, (3.10) 



^ The spin-dependent potentials for quarks of masses mi, 7712 and spin operators Si, S2 read 

Sl • S2 , . , . 1 fx • Sl X ■ S2 



, , 1 rx • Sl X ■ S2 Sl • S2] , - 

V4(r)-\ 5— VbW. 
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is needed including the 2-loop term, in order to have a neghgible perturbative uncer- 
tainty. We calculated the coefficient in the expansion 

7^^{gsF) = -5sf(7o + 7f 5sF + • • •) (3-11) 

by conversion from the MS scheme to the SF scheme. The 2-loop anomalous dimension 

in the MS scheme, 

7^ = - iiiVf)/(32^^) , (3.12) 

is known from [7, 23-25] . 

For the relation between the two schemes, we use the connection of the operators 
in the lattice minimal subtraction scheme ("lat") and the MS scheme of dimensional 
regularization on the one hand and performed a 1-loop computation of the SF renor- 
malization factor in the lattice regularization on the other hand. The latter, new, 
computation is detailed in App. [Al Here we just quote the results and combine them 
to obtain 7^^. 

At 1-loop order, the operator in the lattice minimal subtraction scheme is given by 
O'^M = [1 - gho Hfia)]Ospin (3.13) 

with the bare operator Ospin defined earlier. It is related to the operators in the SF 
scheme and the MS scheme by finite renormalizations, 

Ol^fi) = XSF,lat(5lit(/^))Oipin(/^) , = XMSMt(sfM)OivM ■ (3-14) 

For our purposes it suffices to know the expansions up to order 5^, 

Xa,6(5') = l + x25' + ... , (3.15) 

where 

= 0.3824 (3.16) 



has been computecjl by Flynn and Hill [8] and 



^ 0.3187016(1) - 0.027448(1) Nf (3.17) 



is obtained in App. |Al These are combined to 



^SF,MS = ^SF^iat " xS,iat = "0-0637 " 0.0275 N, . (3.18) 



^Taking into account the discussion in [32], the value of e should be reduced from e = 24.48 in [8] to 
e = 4.53. In the notation of [8] we have = t^[C/(I>c -e + 4-K'^) + CA{Da + Di, - Dc/2)]. We 

thank Jonathan Flynn for clarifying this point. When numerical uncertainties are not written explicitly, 
they are estimated to be at most of order 2 on the last digit. 
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Note that the A'f -dependent part depends on the value 9 chosen for the spatial boundary 
conditions of the quark fields in eq. (j3.8p . The above result refers to 6 = — vr/3, for which 
the expectation values in eq. p.Sp can be shown to be real [33]. 

The last missing ingredient is the relation between the couplings, (at the same 
renormalization scale) 

with the 1-loop coefficient [27, 34] 

X^^ = -^(ci,o + ci,iA^f) , ci,o = 1.25563(4) , ci,i = 0.039863 (3.20) 

for the SF-coupling defined at 6 = ir/b [34]. 

Analogous to [35], where the discussion is carried through for the anomalous di- 
mension of quark masses, we then obtain the desired result 

^SF = + 2bo x'il^ - lox'i^ = -0.00236 - 0.00352 + 0.00023 . (3.21) 

Our computation used the lattice regularization but the result is regularization inde- 
pendent. 



4 Step scaling functions 
4.1 Definition 

The anomalous dimension, which we just obtained in 2-loop approximation, describes 
the change of the operator Ofj^^ under an infinitesimal change in the renormalization 
scale. For numerical, non-perturbative computations one considers finite changes of the 
scale, typically by a factor of two. These define the step scaling function (Tspin(w) via 

0!Uf^) = aspin(5'(l//^)) 0!^,^i2fi) . (4.1) 

It is given by the continuum limit 

CTspiniu) = lim i;spin('U,a/-^^) (4.2) 

a/L^O 



of the lattice step scaling function 



^fp^n(2i) 



(4.3) 



Here g'^{L) denotes the SF coupling as before and the condition m = means that 
the renormalization condition is imposed at vanishing quark mass. In our numerical 
implementation the latter does not play a role since we work in the pure gauge theory. 
It is important to understand whether a given renormalization condition leads to small 
or large a-effects. This determines whether the limit in eq. ()4.2p can be taken by 
an extrapolation from the accessible lattices. A first understanding can be sought in 
perturbation theory. 
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4.2 Lattice artefacts in perturbation theory 

We define the relative lattice artefacts as 

6iu, a/L) = ^^Pin(^.«/-^) -^spin(^) ^ ^^/^^ ^ ^ Q^^2) _ (4 4) 

A term 5o{a/L) is absent since we defined ^^^^^^^(L) such that it is one at tree-level for 
any value of a/L. The 1-loop term may be expanded in A^f, 

6i{a/L) = 6i,o{a/L) + Ni6i,i{a/L) . (4.5) 

With all improvement terms, including the boundary improvement term Ct [27], set to 
their proper perturbative values, the 1-loop cutoff effects turn out to be rather small. 
We show them in Table [T] for the Eichten-Hill action for the static quarks {W{x,0) = 
U{x,0)), the plaquette gauge action and the 0(a)-improved fermion action [20]. 



L/a 


5ifi{a/L) 


5i,i(a/L) 


6 


-0.000236 


0.013742 


8 


-0.000165 


0.005791 


10 


-0.000106 


0.003026 


12 


-0.000072 


0.001876 


14 


-0.000051 


0.001296 


16 


-0.000038 


0.000956 



Table 1: Lattice spacing effects of Egpin in 1-loop perturbation theory, see text. 



4.3 Non-perturbative results for Nf = 

We carried out pure gauge theory simulations to determine Sgpin for different couplings 
u, resolutions 1/12 < a/L < 1/6 and also for different discretizations of the HQET 
action and the operator Tables of the numerical results are found in App. [Bj 

In Fig. [2] we show the a-dependence for a few couplings and choices of discretiza- 
tions. Lattice artefacts are moderate in general and a simple continuum extrapolation 
with an ansatz T, = a + (a/L)'^ p, separately for each value of u seems justifiedQ In cases 
where more than one discretization was simulated at the same value of u, a constrained 
extrapolation (common a but discretization-dependent p) was performed. This simple 
analysis yields the continuum step scaling functions in Table [2l 

A good agreement with perturbation theory is seen at weak couplings, say < 2, 
where, however, the non-perturbative results are not accurate enough to distinguish 
between the 2-loop anomalous dimension and the 1-loop one. On the other hand, for 
couplings 5^ ~ 3 perturbation theory breaks down entirely and the inclusion of the 
2-loop anomalous dimension brings the perturbative curves even further away from the 
MC results. 

^We also checked other extrapolations, see e.g. [11,36]. They give compatible results. 
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0.01 0.02 0.03 



(a/L)^ 

Figure 2: Examples of continuum extrapolations of Egpin for u =1.243, 2.77 and 3.48. Filled symbols 
indicate that was defined as F^^, but with the link variables replaced by HYP2 links. The data at 
the largest coupling has ct at 2-loop precision, otherwise it is set to the 1-loop value. 



u crspin(tt) p(^) action _F)t^ n-loop Ct 



0.8873 


1.0235(68) 


0.02(36) 




0.9944 


1.025(11) 


-0.15(53) 


EH 


1.2430 


1.0302(82) 


0.24(42) 




1.3293 


1.043(12) 


-0.66(55) 


EH 


1.5553 


1.0418(89) 


0.27(43) 




1.8811 


1.075(11) 


-0.44(52) 




2.1000 


1.052(14) 


1.15(69) 




2.4484 


1.089(14) 


-0.31(68) 




2.770 


1.099(11) 


1.60(60) 
-0.01(55) 


HYP2 


3.480 


1.300(34) 


5.2(1.6) 


HYP2 


3.480 


1.308(29) 


3.0(1.4) 


HYP2 


3.480 


1.306(25) 


2.1(1.2) 
3.1(1.3) 


HYP2 



Table 2: Continuum step scaling function CTspin and coefficients of (a/L)^ in the continuum extrapo- 
lations. The standard discretization is the HYP2 action of [18], the clover leaf operator Ffj,^ and the 
1-loop value of Ct. Deviations from this rule are indicated: EH refers to the Eichten-Hill static action 
and the label HYP2 in the column means that the links appearing in F^i, are replaced by HYP2 
links. 
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^spin(u) 



1.3 



1.2 



' ' I ' ' ' ' I ' ' ' ' I 
SF scheme, Nf=0 

- 2-loop 1-loop 

- 3-loop 2-loop 




Figure 3: The step scaling function compared to the perturbative prediction obtained by integration 
of the perturbative RG equations, truncating the (5 function at 3(2)-loop order and the anomalous 
dimension at 2(l)-Ioop order. 



5 Non-perturbative scale dependence for A^f = 

Through an iterative apphcation of the (inverse) step scaling function we now determine 
the non-perturbative scale dependence of the operator in the SF scheme. We first 
represent the numerical data of Table [2] as well as the data for the step scaling function 
of the coupling, o"(n) (Table A. 2 of [11] and Table 4 of [27]) by smooth interpolations 

4 

a{u) = M + ^SjU*+^, (5.1) 

1=0 

o'spin(ii) = l + mou + miu^ + m2U^ (5.2) 

for u < 2.8. Here sq = 21n(2)6o > -si = -Sq + 21n(2)6i , mo = ln(2) 70 are fixed by 
perturbation theory, while the other coefficients are free fit parameters. With these 
parameterizations we then solve the recursions 

uq = f{L^i,x) = 3.48 , a{uk+i) = Uk , ff^(2"''Lmax) = Uk , (5.3) 

-, / ( ^ ^ ^SF(2VL^ax) 

WO = i-, Wk+1 = Wk/Crspm{Uk+l) =^ -T — 7777 7 = Wk ■ 

5'SF(l/-'^maxj 

Errors are propagated through the parameterization and recursion and it is checked 
that changing the number of fit parameters in eq. (|5.ip and eq. ()5.2p does not alter the 
results significantly. Next we apply 



$sf(/u) 



- [ 2bog (/i) J exp <^ - dg 



(g) 70 

/3SF(5) bog 



(5.4) 



$ stands for any matrix element of the operator Ospin, for example the matrix element A2 
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Figure 4: Scale dependence of Cspin in the SF scheme. 



for /i = 

max and the perturbative approximations to the /3-function and anomalous 
dimensions and obtain $sF(l/^max)/^RGl- Finally we form $sF(l/(2Lmax))/^RGl = 

C^spin (Mo)$SF(l/^max)/^'RGI with 

o'spin('iio) from Table [2l For not too small A;, the 
final result is independent of it; the use of perturbation theory is safe in that region. 
Taking k = 6, the 2-loop approximation for the anomalous dimension and the 3- loop 
approximation for the /^-function we arrive at 

$sf(/")/^'rgi = 0.759(17), at = l/L^ax , (5.5) 
^SF(/i)/^RGi = 0.992(29), at /i = l/(2Lmax) . (5.6) 

Figure H] compares the non-perturbative running of Ospin to perturbation theory. 



5.1 The relation of bare and renormalization group invariant operator 

The universal result eq. (15.6P has to be combined with values of Zfp^^{L) at L = 2Ljna 
1.436 ro [37] which depend on the bare coupling and lattice action, to form 



for the respective action. The numerical values of Table [3] are well represented by 



Zfp'i„(2Lmax) = 2.58 + 0.14 (/?- 6) -0.27 (/?- 6)2 EH action, (5.8) 
Zfp'i„(2Lmax) = 2.59 + 0.11 (/3 -6)- 0.34 (/?- 6)2 HYP2 action . (5.9) 

These interpolations may be used in the interval 6.0 < (3 < 6.5 with an error of about 
1%. While this error ought to be taken into account before the continuum extrapolation 
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of the renormalized matrix elements, the uncertainty in eq. (j5.6p apphes additionahy in 
the continuum hmit. 



p 


L/a 


action 


^spin 


6.0219 


8 


EH 


2.585(19) 


6.1628 


10 


EH 


2.602(24) 


6.2885 


12 


EH 


2.602(24) 


6.3992 


14 


EH 


2.589(29) 


6.4956 


16 


EH 


2.593(35) 


6.0219 
6.1628 
6.2885 
6.3992 
6.4956 


8 

10 
12 
14 
16 


HYP2 
HYP2 
HYP2 
HYP2 
HYP2 


2.593(21) 
2.589(26) 
2.589(26) 
2.585(21) 
2.553(22) 



Table 3: Renormalization factor at the matching scale. In all cases F^^, is the standard clover operator. 



6 First applications 

We illustrate the usefulness of our result with two sample applications. 
6.1 Spin splitting 

First we take numbers for the bare A2 which have been reported in the literature. 
Unfortunately they exist only for [3 = 6.0, which corresponds to a ~ 0.1 fm. The more 
recent evaluations are (the light quark has the mass of the strange quark) 



Ref. [5]: a^A^^''^ = 0.0100(19) , (6.1) 
Ref. [6]: a^A^^''*' = 0.0138(15) . (6.2) 

The authors of [5,6] then estimate the mass splitting as 

Ref. [5]: Am^ = 0.28(6) (?) GeV^ , (6.3) 
Ref. [6]: Am^ = 0.36(4) (?) GeV^ , (6.4) 

where the renormalization factor is taken from perturbation theory using a boosted 
coupling [4] and by tadpole improved perturbation theory [6^. With Cspin(Mb/A^) = 
1.15, Mb = 6.76(9) GeV [13, 16] and Z^f^ = 2.6 , a = 1/(2 GeV) we find from eq. dUl) 

Ref. [5] and NP Z^f^: An? = 0.38(7) (?) GeV^ , (6.5) 
Ref. [6] and NP Zf^l: Am^ = 0.53(6) (?)GeV2 . (6.6) 



®Also somewhat different values for the lattice spacing were used by the two groups. 
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In all these estimates the additional uncertainty marked as (?) refers to lattice artefacts, 
namely the fact that a continuum limit has not been taken, and of course to the missing 
dynamical quark determinant. The experimental mass splitting is Am^ = 0.497 GeV^. 

6.2 Renormalization factor for spin-dependent potentials 

In phenomenological applications of the spin-dependent potentials, the standard renor- 
malization scheme is MS, as in [7,23-25]. We apply eq. (j5.4p in the MS scheme with 
A'f = 0, = 238 MeV [11], the 4-loop /3-function and the 2-loop anomalous dimen- 
sion. This yields 



where the cited error bar is half of the change when one uses the 1-loop anomalous 
dimension instead. Remember from Sect.[3]that the square of this renormalization factor 
enters the potentials. As an illustration, for the standard Eichten-Hill action for the 
static quarks, at 50 = 1 (/? = 6.0), we then have Z^{2GeY) « 0.756 x 2.58x0.99 = 1.93 
or a somewhat smaller number at larger fj,. One can compare this at go = 1 to the 
renormalization factors used in [31]. A tree-level tadpole-improved factor is Ztad = 1.684 
and the Huntley-Michael factoi|I| is around Zhm = 1.62. Since a renormalization scheme 
and scale are not specified in these procedures, there is no reason to expect a closer 
agreement. 

7 Conclusions 

We have presented yet another example that the non-perturbative renormalization pro- 
gramme using recursive finite size techniques [9, 11] can be carried out also in difficult 
cases. Four-fermion operators were renormalized successfully in [39,40], and those con- 
taining static quarks in [41,42]. Here we have treated the case of an operator with static 
quarks and gluon fields. 

In several cases quite significant deviations from perturbation theory had already 
been observed at intermediate to low renormalization scales [39,43,44], but the present 
case is the strongest example in that respect (Fig. [3l Fig. [3]). 

Our example in the previous section illustrates that the non-perturbative Z^^^ 
has a rather big effect. Although the quenched estimates for Am^ = mg, — rn^^ are 
in rough agreement with the experimental mass splitting Am^ = 0.497 GeV^ after 
the non-perturbative is used, it now remains to improve the precision of the 

bare matrix element as well as to obtain it at smaller lattice spacings in order to see 
whether the quenched approximation does indeed give a reasonable estimate of the 

^The prescription of [38] yields a small r-dependence, which is not present in the standard renormal- 
ization of local operators. 




(6.7) 
(6.8) 
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spin splitting. Wc also emphasize that there is a remaining uncertainty in the use of 
(continuum) perturbation theory for Cspin(Mb/Aj^). This can be significantly reduced 
by a computation of the associated 3-loop anomalous dimension, but also an entirely 
non-perturbative matching of HQET and QCD is promising [12,13]. 
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A One-loop computation 

Our aim is to compute the expectation value of a general Wilson loop at one-loop order 
in the Schrodingcr functional (SF), bearing in mind that we are finally interested in the 
computation of the expectation value of a clover operator F^i, inserted into a Polyakov 
loop, enabling us to obtain the two-loop anomalous dimension of the chromo-magnetic 
operator in the SF scheme from its value in the lat scheme. Due to the space-time 
locality of such an observable, it will be advantageous to compute the gluon loops in x 
space, while the contribution of tadpoles is proportional to the zero-momentum gluon 
propagator. Quite a lot of notation will be introduced, but the formulae presented here 
are suitable for the automated computation of arbitrary Wilson loops. 

The lattice spacing is set to one in this appendix. Space-time indices run from 1 
to 4, the latter being associated with time. For everything else we reuse the notation 
of [45], referred to as (K) in what follows, except that a twiddle on the color components 
of gluon, ghost and quark fields is dropped. The reader is assumed to be familiar with 
chapters 3, 4 and 5 of (K). Up to one-loop order, the observable O has the perturbative 
expansion (K:4.39) 

{O) = + ai _ (O(i)5(i))o] + 0{gt), (A.l) 

where (. . .)o is the expectation value with respect to the free part S^^^ of the action and 
O^'^) is defined by O = O^^) + O^^^go + O^'^^g^ + .... 

The (constant, Abclian) background field induced by the non-trivial boundary con- 
ditions takes the value V"^(a;). We use the basis of the su{3) Lie algebra 

8 
a=l 
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defined in (K:App. A) as well as the Fourier representation 

<lt{^) = ^E^"'"«4(P,^4). (A.3) 
p 

Iki^) = ^E"'''"«'^'"^'^'^^"'^^^'«fe(P'^4). (A.4) 
p 

The gluon propagator in mixed representation used by (K) has the form 

(9^(P, 2:4) qt{p', yi))o = SbaL^5p+p> D'j^^{p; X4, Va)- (A.5) 

Kronecker symbols such as Sp or (5ju_4 carrying a single index are shorthands for Sp^o 
and (5^^4 respectively. 

A.l Parameterization of the observable 

In order to compute the expectation value of an arbitrary Wilson loop at one-loop order, 
we parameterize the loop by a starting point x^^^^^^^ and an ordered list £ of length £. 
The entries of the list are directions /if, i = !,...,£. These directions take non-zero 
integer values between —4 and +4. An electric plaquette in the (03) plane is thus 
parameterized by £ = (3 4 — 3 — 4). Clearly the loop is closed if and only if each 
integer appears as many times with the + sign as it docs with the — sign (modulo L for 
the space directions). We normally drop the £ in /if since we will be dealing only with 
one path at a time. 

The sequence of points the loop goes through is obtained as follows, 

^(1) = ^(start) ^ ^(i) ^ ^.^ i = 1, 1. (A.6) 

fi = sign(/i)|/i| are unit vectors pointing in the four it directions of the lattice. We 
identify 4 = 0. At tree-level, the expectation value of the Wilson loop is 

e 

W^V] = llv{x^^\^i,). (A.7) 

i=l 

In general, for any 4-vector we introduce negative-index components 

P-^c = -Pix- (A.8) 

Because of the way the path is parameterized, for any link variable we introduce 
negative-index components by imposing 

U{x,tJ,) = U^{x + il,-tJ,), ql{x) = -q\(x + ii). (A.9) 

The Fourier representation is now defined for all /x as follows, 

= 1^ E e^'''(P'-^''^) g^(p, X4) (A.IO) 
p 
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where 

g;i(p,x4) = gf^|(p,x4 - (A.ll) 

and 



(I if = 4 

giiPk+M^'i))/'^ if /X = A; 

_e«{-Pfe+0a(a;4))/2 if ^ = — A: 

-1 if = -4 



= sign(/.) (5|^|_4 + (1 - 5|^|_4)e^(P-+*''(-4))/2) . (A.12) 
With these notations we have 



illiP, X4:)qtip', Va)) = (5fea^^5p+p'-Dp^||^| (p; X4 - (5^+4; ^4 - S^+a) ■ (A.13) 
A.2 Tadpoles: {0^^'^S^^^)o 

The terms considered in this subsection owe their existence to the non-vanishing back- 
ground field. Since the latter is diagonal, the V matrices all commute two-by-two and 
we have 

e 

tr{t^«} = ^ tr{g^^.(x(^)) W^iV]}. (A.14) 
i=i 

The three contributions to S^^^ (coming from the gauge, ghost and quark terms) are 
given by eqs. (K:5.64,5.72) and for the quarks eqs. (K:5.76,5.82). After a short calcula- 
tion one finds 

(tr{w}^)}5W)o = - J2 EE«l;/^4)r;(n4), (A.15) 

ae{3,8}A»=l "4 

with 

e 

al^{u,) = tr {f^WlV]] sign(M,) (5|^,|-4 + (1 " ^,,\-,)e-'^-^^'^''^'' 

i=i 

^;i^,l(o;n4,xi^')-5^^.+4). (A.16) 

and 

(U4). (A. 17) 

^/?,giuon' ^M.ghost' ^/?,quark ^re defined respectively by eqs. (K:5.114), (K:5.115 ) and 
(K:5.116). Note that for the ghost and quark case, the overall minus sign in eq. (|A.15p 
arises because of the loop (ghosts and quarks are anti-commuting) . For the gluons, the 
minus sign is just a convention chosen by (K) and is compensated by a minus sign in 
the definition of T^ gj^on- 
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A.3 Gluon loops: (C'(2))o 

We separately consider two contributions: 

Expanding the exponential of the gluon field to linear order, we obtain 

where we have used the cyclicity of the trace and the shorthand g'--''' = q^.{x^^^). There 
is also a contribution from the quadratic piece of the exponential of the gluon field, 

which is 1/2 of the term j = f in (|A.19P . 

We also need the notation Sj = ^(1 — sign(/ij)) and n = 1 + mod(n— 1, for n > 1. 
Now we can formulate the definition 



WAisM' s;)-l ^^^^ if j + + and s,=0 

W^J s,\j sy) - j ^(j^ ^ ^ otherwise ^ ^ ^ 

that invokes the parallel transporter along the loop from x^^^ to x^^'^: 



1 if J = f 

]^■~,^V{x^^,^x,) ifj<j' (A.22) 



I U'i=,v{x(^\fi,)UiJiV{x(^\f,i) if j>/ 

One then finds 

e £ 8 , 
(tr{VFl^''^})o = -3 E E EE^'^^'''''"''*'''^ ^^'^"^'''''"'^''''^ ''''^ ^ 

(A.23) 

We introduce the propagator completely in x-space, 
A;1,(x; x4, 2/4) = ^ E ^'"'^ e^''"(P'-^''^) e^e.(-p,,4,.) ^a^||^| (p. _ 5^^^. _ ^ 

(A.24) 
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which allows us to write 

£ £ 8 

8 £ 

(tr{^f^)})o = i5]tr{/W,-{F]}5^A;i^.^^, (0;xi^\xi^)). (A.26) 



(A.25) 



a=l i=l 



A. 4 Improvement 



In order to be able to reach the continuum limit with a rate proportional to (1/L)^ our 
observable needs to be improved. Since there are no operators of dimension 6 with the 
same symmetries of Ospim non-vanishing at one-loop order, and with no valence quarks, 
the improvement amounts to compute the additional contributions stemming from the 
volume and boundary counter-terms in the action. The volume term is proportional to 
Csw5 whose tree-level expression, Cgw = 1, enters our observable at one-loop order. It is 
taken into account directly in the quark propagator. The only boundary term needed 
is proportional to the one-loop expression of Ct [27,34]. The corresponding counterterm 
can be expressed as 

(tr{Ty«}5«,)o (A.27) 
with S^^^ given in eq. (5.130) of (K). The explicit expression reads 

(tr{H^f )}4\lb)o = ^cf)[sin(27) +sin(7)]tr{/8H^,-iy]} J^M,--,, (A.28) 

fe=i 

with 

M,-., = j;sign(^,-)(5|M.|-4 + (1 - ^|M.|-4)e-^'^«(^^'')/'J (A.29) 



and 7 = tt/SLT once the "point A" has been chosen. The contribution (|A.28p vanishes 
for the Polyakov loop 'P3i^)\^^—j'/2 'without operator insertion. 

A. 5 Summary 

The expectation value of the Wilson loop at one-loop order is given by 

{tv{W^}) = W^V]+gl [{tv{W^'^% + (tr{iy(^'^)})o (A.30) 

-{tv{WP}S^%-{tr{W^'^}s[ll.) 
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where the one- loop terms are given by eqs. ()A.25IA.26|A.15| ) and the improvement term 
by eq. (IX28|) . 



A. 6 Implementation in MATLAB 

For our perturbative computations we decided to use MATLAB in order to combine 
comfortable programming, robustness of the libraries and acceptable speed for the in- 
volved observables and lattices. 

In presence of a non-vanishing background field, a simple analytical expression for 
the ghost, gluon and quark propagators is not available. They are computed by ex- 
ploiting the recursive techniques presented in [46,47]. The gluon propagator is the most 
time consuming computation. Its Fourier transformed expression ()A.24p is calculated by 
summing only over a reduced set of momenta, which saves a factor of 6 (asymptotically 
on large lattices) in computing time. 

The tadpole loops are observable independent, and they are computed and stored. 
We use the formulae of [48] , where the symmetries of vertices and propagators are fully 
exploited. Then the contributions (jA.lSp are computable with an effort negligible in 
comparison to the loops. 

The improvement counterterm involves only the zero momentum gluon propaga- 
tor and the trace of the product of diagonal matrices; it is computationally cheap in 
comparison to the rest. 

In order to give an idea of the computational cost, for L = 48 the computation 
of all diagrams and improvement counter-terms for the Polyakov loop with insertion 
of the clover leaf operator has been carried out in 2 weeks on a PC, equipped with a 
single processor Intel Pentium 4 with 2.6 GHz. The scaling can be approximated with 
a polynomial in L, and is asymptotically dominated by the highest power, i.e. 

A. 7 Checks 

Of all checks we did to confirm the correctness of our code, we briefly report about two 
of them, which may be of interest in other applications. 

As observed in [49], the expectation value of the gauge action can be evaluated 
economically by taking the logarithmic derivative of the partition function with respect 



The term v arises from the gauge degrees of freedom that are constant in space and live 
on the lower temporal boundary (there is no extra gauge degree of freedom associated 
with the boundary 2:4 = T because global symmetries are not to be gauge-fixed): 




(A.31) 
(A.32) 




with boundary links set to unity 

with non-trivial Abelian boundary field 



(A.33) 
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Figure 5: One- loop connection between the SF scheme and the lat scheme. 



L 


ASF, lat 


(1,1) 

XsF.lat 


4 


0.319032402694607 


-0.071383603501862 


6 


0.319029788544510 


-0.043758342410953 


8 


0.318915972669536 


-0.034594298272846 


10 


0.318837839980287 


-0.031316158100235 


12 


0.318793467352611 


-0.029898520735474 


14 


0.318767160420896 


-0.029160440246041 


16 


0.318750540373758 


-0.028720855529176 


18 


0.318739443194014 


-0.028434643639671 


20 


0.318731691659922 


-0.028236618214183 


22 


0.318726074540605 


-0.028093427502346 


24 


0.318721879237185 


-0.027986322197915 


26 


0.318718665931215 


-0.027904013713480 


28 


0.318716151870336 


-0.027839341836381 


30 


0.318714148841999 


-0.027787573356721 


32 


0.318712527751079 


-0.027745471721788 


34 


0.318711197727617 


-0.027710760078004 


36 


0.318710093332871 


-0.027681797428228 


38 


0.318709166494586 


-0.027657376001968 


40 


0.318708381269461 


-0.027636590425178 


42 


0.318707710343347 


-0.027618750920289 


44 


0.318707132670872 


-0.027603324328397 


46 


0.318706631831595 


-0.027589893181705 


48 


0.318706194851218 


-0.027578126762064 



Table 4: Results for the connection between the SF and the lattice MS schemes. 
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In the first case, all 8 zero momentum gluons at the lower temporal boundary obey 
Dirichlet boundary conditions. They are associated with spatially constant modes and 
are therefore not propagating degrees of freedom. With a non-trivial Abelian back- 
ground field, only two of the gluon fields obey the Dirichlet boundary conditions at the 
lower boundary, and are associated with spatially constant diagonal modes; the others 
are propagating modes. We checked that our program reproduces this result. 

A further successful check, which we do not report in detail, consists in comparing 
the perturbative results for the plaquette and the Polyakov loops, with and without 
insertion of the clover leaf operator, to the corresponding non-perturbative (quenched) 
computations. The latter are performed at small bare couplings, 0.015 < Qq < 0.06, 
setting all needed improvement coefhcients to their tree-level values. In all cases L = 
r = 4. 



A. 8 The Polyakov loop and chromo-magnetic operator 

A. 8.1 Tree-level computation 

As far as the gauge boundary values and the induced background field are concerned, we 
follow [27], and work with the boundary fields defined by eqs. (|3.1|3.2p . The numerator 
and denominator on the r.h.s. of the renormalization condition (13. 8p assume the compact 
form 

3 

L'{i,{V^{x)E^{x)))g^^^ = L2^exp{i;[x4(A™ + (L-X4)<^J} 

(A.34) 

X sin {(j)'^ -(t>m)]- 

3 

(trp3(2;))^^^Q = J] exp{|;[j;4(/>'^ + (L-X4)0J} (A.35) 

m=l 



A. 8. 2 One-loop order 

The one- loop contribution to the lattice step scaling function Sspin(u, 1/L) = 1 + 
E^pjj^(l/L)ii + 0(^2), defined in eq. (g^D, is computed for A^f = and Af = 2 giv- 
ing the results represented in Figure [H There the effect of the 0(a)-improvement is 
evident, especially in the Af = case. In both the unimproved and improved cases the 
continuum limit is consistent with the prediction 7oln(2). 

These results enter the computation of the connection between the SF scheme and 
the lat scheme, as shown in Figure [5l The one-loop connection factor ^ax. obtained 

from the one-loop contribution .^spj^ to the renormalization factor (|3.8p by subtracting 
the logarithmic divergent part 

xi^MiL) = zII^'\l) - 70 ln(L) . (A.36) 
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Figure 6: One- loop contribution to the step scaling function of the chromo- magnetic operator in the 
Ni ^0 (left) and iVf = 2 (right) cases. 

We decompose Xsf lat according to its A'f-dependence 

(1) _ (1,0) AT ^(1.1) /A 07^, 

The quarks are massless, implemented at this order in perturbation theory by mo = 
and the angle 9 = — vr/S. The fuh hsting of results is shown in Tab. [H where all 
numbers are given with 15 digits for readability, although the last two or three may be 
insignificant. 

The continuum limits, including the estimate of the associated uncertainties, are 
performed according to the method described in [47] , with MATLAB routines provided 
by Ulli Wolff. We have verified that the roundoff errors as well as the errors in ct quoted 
in [27, 34] are negligible compared to the systematic uncertainties of the extrapolation. 
The final result of the continuum limit extrapolations is expressed in eq. ()3.17p . 

B Monte Carlo simulations 

In our measurements of observables we fully exploit translational and axis exchange 
invariance. The ensemble of gauge configurations is generated by means of the "hybrid 
over-relaxation" algorithm with lexicographically ordered sweeps (see e.g. [50] for the 
exact implementation). The basic update consists of 1 heat-bath update sweep [51- 
53], followed by A^'qr over-relaxation sweeps [54]. The update is iterated A^up = 2 
times between measurements and the parameter A^oR varies from a minimum of 3, for 
L/a = 6, to a maximum of 10, for L/a = 24. This guarantees to have short integrated 
autocorrelation times for our observables, while the computing time spent for the update 
does not exceed the one required for the measurements. Still there can be very slow 
modes in the system as will be discussed at the end of the appendix. 
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Figure 7: Autocorrelation pattern of Zapin(2L) for two quenched simulations with Lja = 16. Left: 
Normalized autocorrelation function plotted vs. the separation of measurements t. Right: Integrated 
autocorrelation time rint vs. the summation window W. With fi — 7.6547 we have g^{L) = 1.8811 (red 
asterisks), while for /3 = 6.4527 we have g'^{L) = 3.480 ( green circles). The vertical lines correspond to 
the optimal values of W computed according to [57]. 

In the case of the Eichten-Hih action, the gauge hnks building up the Polyakov 
loop, but not the inserted clover leaf operator, are evaluated by a 10-hit multi-hit 
procedure [55], where each hit consists of a heat-bath update of the above type. With 
this variance reduction, the statistical precision is similar to the one of the HYP2 action. 
On the largest lattice {L/a = 24) we could obtain an around 1% precision in Zf^^^{L) 
with 30k measurements at /3 ~ 11 and with 100k measurements at /3 ~ 7.2. However, 
at /? ~ 6.8, where the length of of the Polyakov loop amounts to 0.7 fm, it became very 
costly to reach even a 2% precision. We tried various ways to reduce the variance, in 
particular different versions of multilevel algorithms inspired by [56] , but did not succeed 
in finding significant gains. We then changed the discretization of Ospim replacing the 
links in F^if by HYP2 links. (This is indicated throughout the paper as the discretization 
with HYP2 action and HYP2 operator F^i,.) The resulting increase in precision allowed 
to obtain the last entry in Tab. [5] with 250k measurements. 

As a check that the change of discretization does not introduce unwantedly large 
effects, we also repeated the computation of (t(2.77) this way. Figure [2] nicely confirms 
the expected universality and effects actually turn out to be reduced! We point 
out that with the Eichten-Hill action the cutoff effects can be directly compared to the 
expectations of perturbation theory. For the investigated couplings, g'^ = 0.9944, 1.3293, 
the agreement is very good. 

The large amount of statistically independent measurements needed imposes strong 
limitations to the application of this method to the theory with dynamical fermions. 
The raw simulation results are reported in Table O 

We finally add a remark on error estimates and autocorrelations in our simulations. 
The autocorrelation function F of ^spim defined as in [57], falls very quickly for all our 
mesurements. At all but the largest coupling the integrated autocorrelation time is 
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/3 


L/a 


action 






■^flSn(2-I') 


^spin 


0.8873(5) 


10.7503 


6 






1.3188(46) 


1.3504(50) 


1.0239(52) 


0.8873(10) 


11.0000 


8 






1.3218(22) 


1.3532(51) 


1.0238(42) 


0.8873(30) 


11.3384 


12 






1.3268(30) 


1.3580(69) 


1.0235(57) 


0.9944(7) 


10.0500 


6 


EH 




1.3651(44) 


1.3905(76) 


1.0186(64) 


0.9944(13) 


10.3000 


8 


EH 




1.3514(52) 


1.3924(88) 


1.0303(76) 


0.9944(30) 


10.6086 


12 


EH 




1.3608(53) 


1.384(12) 


1.0171(96) 


1.2430(13) 


8.8997 


6 






1.4336(44) 


1.4839(65) 


1.0351(55) 


-1 C\ A O / -I A\ 

1.2430(14) 


9.1544 


8 






1.4278(30) 


1.4803(56) 


1.0367(45) 


1.2430(35) 


9.5202 


12 






1.4349(27) 


1.474(10) 


1.0275(74) 


1.3293(12) 


8.6129 


6 


EH 




1.4727(57) 


1.5116(77) 


1.0264(66) 


1.3293(21) 


8.8500 


8 


EH 




1.4664(71) 


1.503(12) 


1.0248(95) 


1.3293(60) 


9.1859 


12 


EH 




1.4528(65) 


1.517(13) 


1.0438(99) 


1.5553(15) 


7.9993 


6 






1.5302(61) 


1.6036(48) 


1.0480(53) 


1.5553(24) 


8.2500 


8 






1.5229(46) 


1.5999(78) 


1.0506(60) 


1.5553(70) 


8.5985 


12 






1.5191(27) 


1.579(11) 


1.0394(77) 


1.8811(22) 


7.4082 


6 






1.6193(74) 


1.7208(65) 


1.0627(63) 


1.8811(28) 


7.6547 


8 






1.6100(55) 


1.717(13) 


1.0664(91) 


1.8811(."i8) 


7. !)!)!):> 


12 






l.()()15( 12) 


1.718(1 i) 


1.072()(9i) 


2.1000(39) 


7.1214 


6 






1.6563(87) 


1.792(12) 


1.0819(92) 


C\ '\ f\f\f\ { A 

2.1000(45) 


7.3632 


8 






1.652(11) 


1.773(11) 


1.0732(96) 


O 'X f\f\f\ f C>f\\ 

2.1000(80) 


7.6985 


12 






1.6577(91) 


1.751(17) 


1.056(12) 


2.4484(37) 


6.7807 


6 






1.7618(88) 


1.900(13) 


1.0786(90) 


C\ A A Ci A { A 

2.4484(45) 


7.0197 


8 






1.7371(86) 


1.898(17) 


1.092(11) 


2.4484(80) 


7.3551 


12 






1.7141(89) 


1.855(17) 


1.082(12) 


2.770(7) 


6.5512 


6 






1.8317(75) 


2.085(19) 


1.138(11) 


z.77U(7) 


6.7860 


8 






1.8067(94) 


2.044(17) 


1.131(11) 


2.770(11) 


7.1190 


12 






1.7975(90) 


2.000(23) 


1 110/1 ,i\ 

1.113(14) 


2.770(7) 


6.5512 


6 




HYP2 


1.3659(36) 


1.501(12) 


1.0986(91) 


2.770(7) 


6.7860 


8 




HYP2 


1.3643(24) 


1.505(15) 


1.103(11) 


2.770(11) 


7.1190 


12 




HYP2 


1.3668(63) 


1.490(14) 


1.090(11) 


3.480(8) 


6.2204 


6 




HYP2 


1.4329(53) 


2.070(29) 


1.444(21) 


3.480(14) 


6.4527 


8 




HYP2 


1.4350(89) 


1.975(31) 


1.376(23) 


3.480(39) 


6.7750 


12 




HYP2 


1.4465(61) 


1.937(41) 


1.339(29) 


3.480(8) 


6.257 


6 






1.9864(75) 


2.714(32) 


1.366(17) 


3.180(8) 


(). 176 


8 






1.9 192(75) 


2.()08(:',:!) 


1.:!:58(18) 


3.480(8) 


6.257 


6 




HYP2 


1.4297(47) 


1.979(26) 


1.384(19) 


3.480(8) 


6.476 


8 




HYP2 


1.4262(29) 


1.961(25) 


1.375(17) 


3.480(9) 


6.799 


12 




HYP2 


1.4280(33) 


1.864(37) 


1.305(26) 



Table 5: Raw simulation results. The standard discretization is the HYP2 action of 
[18] and the clover leaf operator Ff^i,. Deviations from this rule are indicated. The 
improvement coeflicient ct is set to its 1-loop value, except for the last five lines, where 
2-loop precision is used. The renormalized coupling is reproduced from [11]. 
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then easily estimated. However for L ~ 0.7 fm, we observe that F shows a long tail 
before approaching zero. This pattern is absent in the smaller volumes, and shows little 
sensitivity to changes of A^'or- Figure [7] (left) plots the F-function obtained from runs 
with (L/a, A^uP; -^or) = (16,2,4), for two different physical volumes. Figure [7] (right) 
shows that for /3 = 7.6547 the integrated autocorrelation time is quite short, i.e. Tint = 
0.86(9), whereas for /3 = 6.4527 the tail mentioned above leads to Tint = 2.15(13). The 
latter translates into an increase of the error by a factor of two, in comparison to the 
case where no correlation is present (i.e. Tint = 0.5). 
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